FixNCut: single-cell genomics through reversible tissue fixation and dissociation.

The use of single-cell technologies for clinical applications requires disconnecting sampling from downstream processing steps. Early sample preservation can further increase robustness and reproducibility by avoiding artifacts introduced during specimen handling. We present FixNCut, a methodology for the reversible fixation of tissue followed by dissociation that overcomes current limitations. We applied FixNCut to human and mouse tissues to demonstrate the preservation of RNA integrity, sequencing library complexity, and cellular composition, while diminishing stress-related artifacts. Besides single-cell RNA sequencing, FixNCut is compatible with multiple single-cell and spatial technologies, making it a versatile tool for robust and flexible study designs.

HAMSAB diet ameliorates dysfunctional signaling in pancreatic islets in autoimmune diabetes.

An altered gut microbiota is associated with type 1 diabetes (T1D), affecting the production of short-chain fatty acids (SCFA) and glucose homeostasis. We previously demonstrated that enhancing serum acetate and butyrate using a dietary supplement (HAMSAB) improved glycemia in non-obese diabetic (NOD) mice and patients with established T1D. The effects of SCFA on immune-infiltrated islet cells remain to be clarified. Here, we performed single-cell RNA sequencing on islet cells from NOD mice fed an HAMSAB or control diet. HAMSAB induced a regulatory gene expression profile in pancreas-infiltrated immune cells. Moreover, HAMSAB maintained the expression of ß-cell functional genes and decreased cellular stress. HAMSAB-fed mice showed preserved pancreatic endocrine cell identity, evaluated by decreased numbers of poly-hormonal cells. Finally, SCFA increased insulin levels in human ß-like cells and improved transplantation outcome in NOD/SCID mice. Our findings support the use of metabolite-based diet as attractive approach to improve glucose control in T1D.

Ligand-dependent hedgehog signaling maintains an undifferentiated, malignant osteosarcoma phenotype.

TP53 and RB1 loss-of-function mutations are common in osteosarcoma. During development, combined loss of TP53 and RB1 function leads to downregulation of autophagy and the aberrant formation of primary cilia, cellular organelles essential for the transmission of canonical Hedgehog (Hh) signaling. Excess cilia formation then leads to hypersensitivity to Hedgehog (Hh) ligand signaling. In mouse and human models, we now show that osteosarcomas with mutations in TP53 and RB1 exhibit enhanced ligand-dependent Hh pathway activation through Smoothened (SMO), a transmembrane signaling molecule required for activation of the canonical Hh pathway. This dependence is mediated by hypersensitivity to Hh ligand and is accompanied by impaired autophagy and increased primary cilia formation and expression of Hh ligand in vivo. Using a conditional genetic mouse model of Trp53 and Rb1 inactivation in osteoblast progenitors, we further show that deletion of Smo converts the highly malignant osteosarcoma phenotype to benign, well differentiated bone tumors. Conversely, conditional overexpression of SHH ligand, or a gain-of-function SMO mutant in committed osteoblast progenitors during development blocks terminal bone differentiation. Finally, we demonstrate that the SMO antagonist sonidegib (LDE225) induces growth arrest and terminal differentiation in vivo in osteosarcomas that express primary cilia and Hh ligand combined with mutations in TP53. These results provide a mechanistic framework for aberrant Hh signaling in osteosarcoma based on defining mutations in the tumor suppressor, TP53.

Systematic benchmarking of single-cell ATAC-sequencing protocols.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Fancm has dual roles in the limiting of meiotic crossovers and germ cell maintenance in mammals.

Meiotic crossovers are required for accurate chromosome segregation and producing new allelic combinations. Meiotic crossover numbers are tightly regulated within a narrow range, despite an excess of initiating DNA double-strand breaks. Here, we reveal the tumor suppressor FANCM as a meiotic anti-crossover factor in mammals. We use unique large-scale crossover analyses with both single-gamete sequencing and pedigree-based bulk-sequencing datasets to identify a genome-wide increase in crossover frequencies in Fancm-deficient mice. Gametogenesis is heavily perturbed in Fancm loss-of-function mice, which is consistent with the reproductive defects reported in humans with biallelic FANCM mutations. A portion of the gametogenesis defects can be attributed to the cGAS-STING pathway after birth. Despite the gametogenesis phenotypes in Fancm mutants, both sexes are capable of producing offspring. We propose that the anti-crossover function and role in gametogenesis of Fancm are separable and will inform diagnostic pathways for human genomic instability disorders.

Inhibition of mutant IDH1 promotes cycling of acute myeloid leukemia stem cells.

Approximately 20% of acute myeloid leukemia (AML) patients carry mutations in IDH1 or IDH2 that result in over-production of the oncometabolite D-2-hydroxyglutarate (2-HG). Small molecule inhibitors that block 2-HG synthesis can induce complete morphological remission; however, almost all patients eventually acquire drug resistance and relapse. Using a multi-allelic mouse model of IDH1-mutant AML, we demonstrate that the clinical IDH1 inhibitor AG-120 (ivosidenib) exerts cell-type-dependent effects on leukemic cells, promoting delayed disease regression. Although single-agent AG-120 treatment does not fully eradicate the disease, it increases cycling of rare leukemia stem cells and triggers transcriptional upregulation of the pyrimidine salvage pathway. Accordingly, AG-120 sensitizes IDH1-mutant AML to azacitidine, with the combination of AG-120 and azacitidine showing vastly improved efficacy in vivo. Our data highlight the impact of non-genetic heterogeneity on treatment response and provide a mechanistic rationale for the observed combinatorial effect of AG-120 and azacitidine in patients.

Epigenetic Activation of Plasmacytoid DCs Drives IFNAR-Dependent Therapeutic Differentiation of AML.

Pharmacologic inhibition of epigenetic enzymes can have therapeutic benefit against hematologic malignancies. In addition to affecting tumor cell growth and proliferation, these epigenetic agents may induce antitumor immunity. Here, we discovered a novel immunoregulatory mechanism through inhibition of histone deacetylases (HDAC). In models of acute myeloid leukemia (AML), leukemia cell differentiation and therapeutic benefit mediated by the HDAC inhibitor (HDACi) panobinostat required activation of the type I interferon (IFN) pathway. Plasmacytoid dendritic cells (pDC) produced type I IFN after panobinostat treatment, through transcriptional activation of IFN genes concomitant with increased H3K27 acetylation at these loci. Depletion of pDCs abrogated panobinostat-mediated induction of type I IFN signaling in leukemia cells and impaired therapeutic efficacy, whereas combined treatment with panobinostat and IFNα improved outcomes in preclinical models. These discoveries offer a new therapeutic approach for AML and demonstrate that epigenetic rewiring of pDCs enhances antitumor immunity, opening the possibility of exploiting this approach for immunotherapies. SIGNIFICANCE: We demonstrate that HDACis induce terminal differentiation of AML through epigenetic remodeling of pDCs, resulting in production of type I IFN that is important for the therapeutic effects of HDACis. The study demonstrates the important functional interplay between the immune system and leukemias in response to HDAC inhibition. This article is highlighted in the In This Issue feature, p. 1397.

Same-Cell Co-Occurrence of RAS Hotspot and BRAF V600E Mutations in Treatment-Naive Colorectal Cancer.

PURPOSE: Mitogen-activated protein kinase pathway-activating mutations occur in the majority of colorectal cancer (CRC) cases and show mutual exclusivity. We identified 47 epidermal growth factor receptor/BRAF inhibitor-naive CRC patients with dual RAS hotspot/BRAF V600E mutations (CRC-DD) from a cohort of 4,561 CRC patients with clinical next-generation sequencing results. We aimed to define the molecular phenotypes of the CRC-DD and to test if the dual RAS hotspot/BRAF V600E mutations coexist within the same cell. MATERIALS AND METHODS: We developed a single-cell genotyping method with a mutation detection rate of 96.3% and a genotype prediction accuracy of 92.1%. Mutations in the CRC-DD cohort were analyzed for clonality, allelic imbalance, copy number, and overall survival. RESULTS: Application of single-cell genotyping to four CRC-DD revealed the co-occurrence of both mutations in the following percentages of cells per case: NRAS G13D/KRAS G12C, 95%; KRAS G12D/NRAS G12V, 48%; BRAF V600E/KRAS G12D, 44%; and KRAS G12D/NRAS G13V, 14%, respectively. Allelic imbalance favoring the oncogenic allele was less frequent in CRC-DD (24 of 76, 31.5%, somatic mutations) compared with a curated cohort of CRC with a single-driver mutation (CRC-SD; 119 of 232 mutations, 51.3%; P = .013). Microsatellite instability-high status was enriched in CRC-DD compared with CRC-SD (23% v 11.4%, P = .028). Of the seven CRC-DD cases with multiregional sequencing, five retained both driver mutations throughout all sequenced tumor sites. Both CRC-DD cases with discordant multiregional sequencing were microsatellite instability-high. CONCLUSION: Our findings indicate that dual-driver mutations occur in a rare subset of CRC, often within the same tumor cells and across multiple tumor sites. Their presence and a lower rate of allelic imbalance may be related to dose-dependent signaling within the mitogen-activated protein kinase pathway.

Non-genetic determinants of malignant clonal fitness at single-cell resolution.

All cancers emerge after a period of clonal selection and subsequent clonal expansion. Although the evolutionary principles imparted by genetic intratumour heterogeneity are becoming increasingly clear1, little is known about the non-genetic mechanisms that contribute to intratumour heterogeneity and malignant clonal fitness2. Here, using single-cell profiling and lineage tracing (SPLINTR)-an expressed barcoding strategy-we trace isogenic clones in three clinically relevant mouse models of acute myeloid leukaemia. We find that malignant clonal dominance is a cell-intrinsic and heritable property that is facilitated by the repression of antigen presentation and increased expression of the secretory leukocyte peptidase inhibitor gene (Slpi), which we genetically validate as a regulator of acute myeloid leukaemia. Increased transcriptional heterogeneity is a feature that enables clonal fitness in diverse tissues and immune microenvironments and in the context of clonal competition between genetically distinct clones. Similar to haematopoietic stem cells3, leukaemia stem cells (LSCs) display heritable clone-intrinsic properties of high, and low clonal output that contribute to the overall tumour mass. We demonstrate that LSC clonal output dictates sensitivity to chemotherapy and, although high- and low-output clones adapt differently to therapeutic pressure, they coordinately emerge from minimal residual disease with increased expression of the LSC program. Together, these data provide fundamental insights into the non-genetic transcriptional processes that underpin malignant clonal fitness and may inform future therapeutic strategies.

Antigen-driven EGR2 expression is required for exhausted CD8+ T cell stability and maintenance.

Chronic stimulation of CD8+ T cells triggers exhaustion, a distinct differentiation state with diminished effector function. Exhausted cells exist in multiple differentiation states, from stem-like progenitors that are the key mediators of the response to checkpoint blockade, through to terminally exhausted cells. Due to its clinical relevance, there is substantial interest in defining the pathways that control differentiation and maintenance of these subsets. Here, we show that chronic antigen induces the anergy-associated transcription factor EGR2 selectively within progenitor exhausted cells in both chronic LCMV and tumours. EGR2 enables terminal exhaustion and stabilizes the exhausted transcriptional state by both direct EGR2-dependent control of key exhaustion-associated genes, and indirect maintenance of the exhausted epigenetic state. We show that EGR2 is a regulator of exhaustion that epigenetically and transcriptionally maintains the differentiation competency of progenitor exhausted cells.

CDK4/6 Inhibition Promotes Antitumor Immunity through the Induction of T-cell Memory.

Pharmacologic inhibitors of cyclin-dependent kinases 4 and 6 (CDK4/6) are an approved treatment for hormone receptor-positive breast cancer and are currently under evaluation across hundreds of clinical trials for other cancer types. The clinical success of these inhibitors is largely attributed to well-defined tumor-intrinsic cytostatic mechanisms, whereas their emerging role as immunomodulatory agents is less understood. Using integrated epigenomic, transcriptomic, and proteomic analyses, we demonstrated a novel action of CDK4/6 inhibitors in promoting the phenotypic and functional acquisition of immunologic T-cell memory. Short-term priming with a CDK4/6 inhibitor promoted long-term endogenous antitumor T-cell immunity in mice, enhanced the persistence and therapeutic efficacy of chimeric antigen receptor T cells, and induced a retinoblastoma-dependent T-cell phenotype supportive of favorable responses to immune checkpoint blockade in patients with melanoma. Together, these mechanistic insights significantly broaden the prospective utility of CDK4/6 inhibitors as clinical tools to boost antitumor T-cell immunity. SIGNIFICANCE: Immunologic memory is critical for sustained antitumor immunity. Our discovery that CDK4/6 inhibitors drive T-cell memory fate commitment sheds new light on their clinical activity, which is essential for the design of clinical trial protocols incorporating these agents, particularly in combination with immunotherapy, for the treatment of cancer.This article is highlighted in the In This Issue feature, p. 2355.

γδ T Cells in Merkel Cell Carcinomas Have a Proinflammatory Profile Prognostic of Patient Survival.

Merkel cell carcinomas (MCC) are immunogenic skin cancers associated with viral infection or UV mutagenesis. To study T-cell infiltrates in MCC, we analyzed 58 MCC lesions from 39 patients using multiplex-IHC/immunofluorescence (m-IHC/IF). CD4+ or CD8+ T cells comprised the majority of infiltrating T lymphocytes in most tumors. However, almost half of the tumors harbored prominent CD4/CD8 double-negative (DN) T-cell infiltrates (>20% DN T cells), and in 12% of cases, DN T cells represented the majority of T cells. Flow cytometric analysis of single-cell suspensions from fresh tumors identified DN T cells as predominantly Vδ2- γδ T cells. In the context of γδ T-cell inflammation, these cells expressed PD-1 and LAG3, which is consistent with a suppressed or exhausted phenotype, and CD103, which indicates tissue residency. Furthermore, single-cell RNA sequencing (scRNA-seq) identified a transcriptional profile of γδ T cells suggestive of proinflammatory potential. T-cell receptor (TCR) analysis confirmed clonal expansion of Vδ1 and Vδ3 clonotypes, and functional studies using cloned γδ TCRs demonstrated restriction of these for CD1c and MR1 antigen-presenting molecules. On the basis of a 13-gene γδ T-cell signature derived from scRNA-seq analysis, gene-set enrichment on bulk RNA-seq data showed a positive correlation between enrichment scores and DN T-cell infiltrates. An improved disease-specific survival was evident for patients with high enrichment scores, and complete responses to anti-PD-1/PD-L1 treatment were observed in three of four cases with high enrichment scores. Thus, γδ T-cell infiltration may serve as a prognostic biomarker and should be explored for therapeutic interventions.See related Spotlight on p. 600.

Prevalence and potential biological role of TERT amplifications in ALK translocated adenocarcinoma of the lung.

AIMS: The advent of specific ALK-targeting drugs has radically changed the outcome of patients with ALK translocated non-small-cell lung cancer (NSCLC). However, emerging resistance to treatment with ALK inhibitors in these patients remains a major concern. In previous studies, we analysed two ALK+ patient cohorts (TP53 wild-type/TP53 mutated) in terms of copy number alterations. All patients belonging to the TP53 wild-type group had mainly genetically stable genomes, with one exception showing chromosomal instability and amplifications of several gene loci, including TERT. Here, we aimed to determine the prevalence of TERT amplifications in these ALK+ lung cancer patients by analysing an independent cohort of 109 ALK translocated cases. We further analysed the copy numbers of numerous cancer-relevant genes and other genetic aberrations. METHODS AND RESULTS: The prevalence of TERT amplifications was determined by means of FISH analyses. Copy numbers of 87 cancer-relevant genes were determined by NanoString nCounter® technology, FoundationOne® and lung-specific NGS panels in some of these TERT-amplified samples, and clinical data on patients with TERT-amplified tumours were collected. Our data revealed that five (4.6%) of all 109 analysed ALK+ patients harboured amplification of TERT and that these patients had genetically unstable genomes. CONCLUSIONS: Our preliminary study shows that ALK+ adenocarcinomas should be evaluated in the context of their genomic background in order to more clearly understand and predict patients' individual course of disease.

Combined BRAF, MEK, and CDK4/6 Inhibition Depletes Intratumoral Immune-Potentiating Myeloid Populations in Melanoma.

Combined inhibition of BRAF, MEK, and CDK4/6 is currently under evaluation in clinical trials for patients with melanoma harboring a BRAFV600 mutation. While this triple therapy has potent tumor-intrinsic effects, the impact of this combination on antitumor immunity remains unexplored. Here, using a syngeneic BrafV600ECdkn2a-/-Pten-/- melanoma model, we demonstrated that triple therapy promoted durable tumor control through tumor-intrinsic mechanisms and promoted immunogenic cell death and T-cell infiltration. Despite this, tumors treated with triple therapy were unresponsive to immune checkpoint blockade (ICB). Flow cytometric and single-cell RNA sequencing analyses of tumor-infiltrating immune populations revealed that triple therapy markedly depleted proinflammatory macrophages and cross-priming CD103+ dendritic cells, the absence of which correlated with poor overall survival and clinical responses to ICB in patients with melanoma. Indeed, immune populations isolated from tumors of mice treated with triple therapy failed to stimulate T-cell responses ex vivo While combined BRAF, MEK, and CDK4/6 inhibition demonstrates favorable tumor-intrinsic activity, these data suggest that collateral effects on tumor-infiltrating myeloid populations may impact antitumor immunity. These findings have important implications for the design of combination strategies and clinical trials that incorporate BRAF, MEK, and CDK4/6 inhibition with immunotherapy for the treatment of patients with melanoma.

From whole-mount to single-cell spatial assessment of gene expression in 3D.

Unravelling spatio-temporal patterns of gene expression is crucial to understanding core biological principles from embryogenesis to disease. Here we review emerging technologies, providing automated, high-throughput, spatially resolved quantitative gene expression data. Novel techniques expand on current benchmark protocols, expediting their incorporation into ongoing research. These approaches digitally reconstruct patterns of embryonic expression in three dimensions, and have successfully identified novel domains of expression, cell types, and tissue features. Such technologies pave the way for unbiased and exhaustive recapitulation of gene expression levels in spatial and quantitative terms, promoting understanding of the molecular origin of developmental defects, and improving medical diagnostics.

Reprogramming roadmap reveals route to human induced trophoblast stem cells.

The reprogramming of human somatic cells to primed or naive induced pluripotent stem cells recapitulates the stages of early embryonic development1-6. The molecular mechanism that underpins these reprogramming processes remains largely unexplored, which impedes our understanding and limits rational improvements to reprogramming protocols. Here, to address these issues, we reconstruct molecular reprogramming trajectories of human dermal fibroblasts using single-cell transcriptomics. This revealed that reprogramming into primed and naive pluripotency follows diverging and distinct trajectories. Moreover, genome-wide analyses of accessible chromatin showed key changes in the regulatory elements of core pluripotency genes, and orchestrated global changes in chromatin accessibility over time. Integrated analysis of these datasets revealed a role for transcription factors associated with the trophectoderm lineage, and the existence of a subpopulation of cells that enter a trophectoderm-like state during reprogramming. Furthermore, this trophectoderm-like state could be captured, which enabled the derivation of induced trophoblast stem cells. Induced trophoblast stem cells are molecularly and functionally similar to trophoblast stem cells derived from human blastocysts or first-trimester placentas7. Our results provide a high-resolution roadmap for the transcription-factor-mediated reprogramming of human somatic cells, indicate a role for the trophectoderm-lineage-specific regulatory program during this process, and facilitate the direct reprogramming of somatic cells into induced trophoblast stem cells.

A P53-Independent DNA Damage Response Suppresses Oncogenic Proliferation and Genome Instability.

The Mre11-Rad50-Nbs1 complex is a DNA double-strand break sensor that mediates a tumor-suppressive DNA damage response (DDR) in cells undergoing oncogenic stress, yet the mechanisms underlying this effect are poorly understood. Using a genetically inducible primary mammary epithelial cell model, we demonstrate that Mre11 suppresses proliferation and DNA damage induced by diverse oncogenic drivers through a p53-independent mechanism. Breast tumorigenesis models engineered to express a hypomorphic Mre11 allele exhibit increased levels of oncogene-induced DNA damage, R-loop accumulation, and chromosomal instability with a characteristic copy number loss phenotype. Mre11 complex dysfunction is identified in a subset of human triple-negative breast cancers and is associated with increased sensitivity to DNA-damaging therapy and inhibitors of ataxia telangiectasia and Rad3 related (ATR) and poly (ADP-ribose) polymerase (PARP). Thus, deficiencies in the Mre11-dependent DDR drive proliferation and genome instability patterns in p53-deficient breast cancers and represent an opportunity for therapeutic exploitation.

Renal epithelial cells retain primary cilia during human acute renal allograft rejection injury.

OBJECTIVES: Primary cilia are sensory organelles which co-ordinate several developmental/repair pathways including hedgehog signalling. Studies of human renal allografts suffering acute tubular necrosis have shown that length of primary cilia borne by epithelial cells doubles throughout the nephron and collecting duct, and then normalises as renal function returns. Conversely the loss of primary cilia has been reported in chronic allograft rejection and linked to defective hedgehog signalling. We investigated the fate of primary cilia in renal allografts suffering acute rejection. RESULTS: Here we observed that in renal allografts undergoing acute rejection, primary cilia were retained, with their length increasing 1 week after transplantation and remaining elevated. We used a mouse model of acute renal injury to demonstrate that elongated renal primary cilia in the injured renal tubule show evidence of smoothened accumulation, a biomarker for activation of hedgehog signalling. We conclude that primary cilium-mediated activation of hedgehog signalling is still possible during the acute phase of renal allograft rejection.

Targeting enhancer switching overcomes non-genetic drug resistance in acute myeloid leukaemia.

Non-genetic drug resistance is increasingly recognised in various cancers. Molecular insights into this process are lacking and it is unknown whether stable non-genetic resistance can be overcome. Using single cell RNA-sequencing of paired drug naïve and resistant AML patient samples and cellular barcoding in a unique mouse model of non-genetic resistance, here we demonstrate that transcriptional plasticity drives stable epigenetic resistance. With a CRISPR-Cas9 screen we identify regulators of enhancer function as important modulators of the resistant cell state. We show that inhibition of Lsd1 (Kdm1a) is able to overcome stable epigenetic resistance by facilitating the binding of the pioneer factor, Pu.1 and cofactor, Irf8, to nucleate new enhancers that regulate the expression of key survival genes. This enhancer switching results in the re-distribution of transcriptional co-activators, including Brd4, and provides the opportunity to disable their activity and overcome epigenetic resistance. Together these findings highlight key principles to help counteract non-genetic drug resistance.